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ABSTRACT 

Significant gravitational wave emission is expected from gamma-ray bursts 
arising from compact stellar mergers, and possibly also from bursts associated 
with fast-rotating massive stellar core collapses. These models have in common a 
high angular rotation rate, and observations provide evidence for jet collimation 
of the photon emission, with properties depending on the polar angle. Here 
we consider the gravitational wave emission and its polarization as a function 
of angle which is expected from such sources. We discuss possible correlations 
between the burst photon luminosity, or the delay between gravitational wave 
bursts and X-ray flashes, and the polarization degree of the gravitational waves. 

Subject headings: gravitational waves - polarization - binaries: close - black hole 
physics - stars:neutron - gamma rays: bursts 

1. Introduction 

The observations of gamma-ray burst (GRB) afterglows at energies ranging from X-rays 
to radio have lead to an increased understanding of the possible geometry of the inferred 
relativistic outflow or ejecta (see van Paradijs, Kouveliotou Sz Wijers 2000 for a review; also 
Rhoads 1999; Sari, Piran Sz Halpern 1999). At least for the class of long bursts (durations 
> 10 s) there is now significant evidence that the ejecta is collimated in a jet. 

The evidence for jets is based on a change of the GRB light curve time-dependence power 
law index, and there appears to be a correlation in the sense that bursts with the largest 
gamma-ray fluences have the narrowest implied opening angles. Frail et al. (2001), Piran et 
al. (2001) and Panaitescu Sz Kumar (2002) reported that the total gamma-ray energy release, 
after correcting for the collimation and distance determined by the afterglow observations, 
are narrowly clustered around 5 x 10 50 erg, and suggested that the broad distribution in 
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fluence and luminosity for GRBs is largely the result of a wide variation of the opening 
angles. This result can be obtained assuming that the jet emission inside some angle 8j is 
approximately uniform and then drops off abruptly beyond that. However, Rossi, Lazzati 
& Rees (2002) and Zhang & Meszaros (2002) showed that the observations can also be 
interpreted in terms of an alternative model where the jet, rather than having a uniform 
profile out to some definite cone angle, has a universal angle-dependent beam pattern with 
a luminosity per unit solid angle which is maximal along the axis, and drops off gradually 
away from the axis. In this model, it is the difference of viewing angles 9 which causes a 
wide variation in the apparent luminosity of GRBs, L oc 9~ 2 . 

Gravitational waves (GW) are also expected from some types of GRB (e.g. Kobayashi 
& Meszaros 2002). The GW emissivity and its polarization are angle-dependent, and may 
in principle be measurable, depending on the signal to noise ratio and on the detector align- 
ments. In this Letter we discuss the prospects of exploiting measurements of the gravitational 
polarization degree to obtain information on the nature and orientation of the GRB system, 
providing additional constraints on the luminosity of the bursts, and we mention possible 
implications for the interpretation of dark GRB and X-ray flashes. 



2. Gravitational Waves from GRB Progenitors 

Most GRB models involve a stellar system resulting in a rotating black hole and a 
massive (~ 0.1 — 1M Q ) disk of accreting matter around it. Such massive disks can form from 
the fall-back of debris during the formation of the central compact object, which ultimately 
is likely to lead to a black hole. Several scenarios could lead to a black hole-accretion disk 
system (e.g. Meszaros, 2002). This includes double neutron star binary mergers (Eichler et 
al 1989; Ruffert et al. 1997), neutron star - black hole mergers (Paczynski 1991; Janka et 
al. 1999) and massive stellar collapses (collapsars) (Woosley 1993; MacFadyen & Woosley 
1999). 

Coalescing compact binaries are one of the most promising GW sources detectable by 
the laser interferometer gravitational wave observatory (LIGO), and other detectors such 
as VIRGO, GEO600 and TAMA300. The binary coalescence process can be divided into 
three phases: inspiral, merger and ringdown (e.g., Flanagan & Hughes 1998). (1) During the 
inspiral phase, the gravitational radiation reaction time scale is much longer than the orbital 
period. As the binary losses energy by gravitational radiation, the masses gradually spiral 
in toward each other. (2) The merger begins when the orbital evolution is so rapid that 
adiabatic evolution is no longer a good approximation. The two masses then violently merge 
to form a black hole on a timescale comparable to a small multiple of the dynamical time. A 
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bar-mode instability develops in the early stages of the merger, releasing a fraction of their 
rest mass energy in GWs. However, a significant fraction of the stellar material retains too 
much angular momentum to cross the black hole horizon promptly. This creates a temporary 
disk of debris material around the black hole, whose accretion can power a GRB jet. (3) The 
black hole is initially deformed, and in a ringdown phase radiates away the energy associated 
with these deformations as GWs, until it settles down into a Kerr geometry. 

Massive rotating stellar collapses and their associated GW emission have been calculated 
numerically (e.g. Fryer et al. 1999; McFadyen & Woosley 1999; Rampp et al 1998). These 
numerical estimates are not conclusive, as a number of effects (including general relativity, 
secular evolution, etc, see Rampp et al 1998) have been neglected, but they suggest that 
GW emission from massive collapses may be much less important than from compact binary 
mergers. On the other hand, recent semi- analytical estimates (van Putten 2001; Fryer, 
Holz & Hughes 2002; Davies et al 2002) have indicated that massive collapses might emit 
significantly stronger GW signals than expected from the previous numerical estimates. 
Collapsars, i.e. massive stellar collapses leading to a GRB, require a high core rotation rate 
to form a centrifugally supported disk around a central, possibly spinning black hole. A high 
rotation rate, however, is conducive to the development of bar or fragmentation instabilities 
in the collapsing core or/and in the massive disk (Nakamura & Fukugita 1989; Bonnell & 
Pringle 1995; van Putten 2002). The asymmetrically infalling matter perturbs the black hole 
geometry, which leads to ringdown gravitational radiation. 



3. Polarization of Gravitational Waves 

We have seen that various mechanisms, including binary, bar and fragmentation insta- 
bilities and oscillations of black holes, could occur in GRB progenitor systems, resulting in 
GW emission. It is known that the GWs generated by a binary, bar or fragmentation insta- 
bility are dominated by modes with spherical harmonic index I = m = 2 mode (e.g. Misner, 
Thorne & Wheeler 1978). The most slowly damped quasi-normal mode of a Kerr black hole 
also has indices I — m — 2 (Detweiler 1980). Since this mode may be preferentially excited in 
the presence of binary masses or fragmentation of a massive disk, the ringdown gravitational 
radiation is also dominated by I = m = 2 mode. 

The polarization components of / = m = 2 mode depend on the inclination angle a (e.g. 
Thorne 1980) as 

h+ oc (1 + cos 2 a), /i x oc2cosa. (1) 

GRB progenitors emit GWs more strongly along the polar axis than in the equatorial plane, 
the latter being the orbital plane of the binary, the disk fragments or the bar. Since the 
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GRB jets are launched along the polar (angular momentum) axis (i.e. a = 9), the GRB 
progenitors are stronger GW sources than the average non-bursting merger or collapse, as 
pointed out by Kochanek and Piran (1993) in the context of a neutron star binary model. 
When we observe a GRB from a viewing angle respect to the polar axis, the GW signal is 
stronger by a factor of (5/32) 1//2 [(l + cos 9) 4 + (1 — cos 0) 4 ] 1 / 2 than the value averaged over all 
possible viewing angles. Since for the angle dependent jet model (Rossi et al. 2002; Zhang 
& Meszaros 2002) the difference of the viewing angles causes a variation of the apparent 
luminosity of the GRB, we expect a correlation between the apparent luminosity of GRBs 
and the amplitude of the associated GWs, even after distance corrections. Also in the case 
when the jet profile is uniform inside an opening angle 9j, the typical viewing angle is ~ 2^/3, 
and the correlation is still expected. 

The GWs will be circular polarized when viewed along the polar axis (9 = 0), while the 
+ polarization dominates when viewed along the equatorial plane (9 = ir/2). Therefore, the 
polarization of the GWs is also correlated with the luminosity of the GRBs. We define a po- 
larization tensor in analogy with the electromagnetic approach (e.g. Landau & Lifshitz 1975). 
A plane wave at the position of a detector can be written in the form h + = Re {A + e~ tu)t } and 
h x = Re {A x e~ tujt }. If the plane wave is monochromatic, the complex amplitudes A + and 
A x will be constants, whereas if the wave contains frequencies in a small interval Au; we take 
uj to be some average frequency inside this range. A + and A x are slowly varying functions 
of time, compared to the wave oscillation period. The polarization tensor p a b(ci,b = +, x) 
is defined by using the time averaged values of the amplitudes {A a AD and described by the 
Stokes parameters £i , £ 2 and £ 3 : 

(M! - l ( 1 + ^ (o) 

pab ~ (L4 + | 2 + |A X | 2 > 2U+« 1-6 )' 1 ] 

For the / = m = 2 mode, the parameters are 

8cosfl(l + cos 2 fl) 
U ' ^ 2 (l + cos£) 4 + (l-cos#) 4 ' { ) 

= 2(l + cosfl) 2 (l-cosfl) 2 

?3 ~ (l + cos#) 4 +(l-cosfl) 4 " 

The degree of circular polarization and of linear polarization are given by £2 and \/£ 2 + £f = 
£ 3 = P, respectively. Thus, GRB progenitors whose polar axes are directed towards the 
earth would produce brighter GRBs, and circularly polarized GWs. As the viewing angle 
becomes larger relative to the polar axis, the luminosity of the GRB decreases oc 9~ 2 , while 
the degree of linear polarization increases as P ~ 10~ 2 (#/30 degree) 4 . 
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4. 



Polarization Measurement 



The response of an interferometer (interferometer 1) to the gravitational radiation is 
given by a linear combination of two polarization components mi = F + ^h + + F x ^h x where 
the antenna patterns F +1 and F x l depend on the orientation of the interferometer with 
respect to the GW source (e.g. Finn & Chernoff 1993). We assume that the interferometer 
arms are of the same length and that they meet at right angles. Define a right-handed 
coordinate system with one interferometer arm along the x-axis and the other along the 
y-axis. For simplicity, we assume that the source is in the direction of the z-axis. Since we 
can determine the position of a source in the sky by using observations of the GRB and 
afterglow, it is possible to generalize the following discussion to the case of a source with 
an arbitrary sky position. The angular momentum vector of the source (direction of the 
GRB jet) may be oriented in an arbitrary direction. We assume that the projection of the 
angular momentum vector to the x — y plane makes an angle ( with the x-axis. With these 
conventions, the antenna patterns are given by F +1 = cos2£ and F x l = sin 2(. 

To determine the polarization of GWs, one needs a network consisting of at least two 
interferometers which have different (i.e. non-parallel) arm orientations. Consider an identi- 
cal interferometer (interferometer 2) at the same location as interferometer 1. (If we set the 
arrival time of the GRB signal as the origin of time at each interferometer, we can correct 
for the actual physical separation so that the two interferometers can be always considered 
to be at the same location). We assume that the interferometer 2 is rotated by an angle 
— 7r/4 around the z-axis with respect to the interferometer 1. The response of the interfer- 
ometer 2 is m 2 = F +>2 h + + F X)2 h x , and the antenna patters are given by F +>2 = — sin 2^, 
and F x 2 = cos 2^. 

The detection of the polarization of GWs requires observations with a high signal-to- 
noise ratio (SNR) p. A detection is likelier in an optimal case where the wave forms of the 
polarized components, / + and f x , are known. We define the noise- weighted inner product 

as 



where Sj(/) denotes the Fourier transform of the outputs of the two interferometers Si(t) = 
mi(t) + rii(t) (i — 1, 2), n^t) is the noise of the interferometers, S h (f) is the one-sided noise 
power spectral density and f a is the Fourier transform of f a (a = +, x). We normalize the 
functions f a as (f a ,fb) — $ab, hence the dispersion of the noise is unity o 2 = (m, f a ) 2 = 1 
where the overline represents ensemble average. The linear polarization degree P = £3 is 




(5) 



P 



(mi, f+f + (m 2 , f+f - K /x) 2 ~ K fx} 2 
(mi, f+} 2 + (m 2 , f+) 2 + (mi, f x ) 2 + (m 2 , f x ) 2 



(6) 
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Since the outputs of the interferometers Si are always superpositions of signals m, and 
noises n,, the inner products (m*, /„) themselves are not observable. We estimate (m^ f a ) 
by (m„ / ) ~ (sj, f a ). Substituting this into eq. (6), we can get an estimate P' of the polar- 
ization degree. This approximation introduces an error A in the estimate of the polarization 
degree P, given by 

g-P = l + (A-B)/P 

P l+A + P ' V ^ 

A = 2a+(Ai+cos2C - A 2+ sin2C) + A 2 + + \ 2 2+ , (8) 
B = 2a x (A lx sin2C + A 2x cos2C) + A? x + A 2 X , (9) 

where a+ = >/2(l + cos 2 0)[(1 + cosfl) 4 + (1 - cosO) 4 }- 1 / 2 , a x = 2 v / 2cos^[(l + cosfl) 4 + (1 - 
cos#) 4 ]~ 1/2 , \ ia = (m, f a ) lp and the SNR p is p 2 = (h+, f+f + (h x J + ) 2 . The fluctuations 
\ ia are of order of p~ x . Thus, a SNR of p > P -1 is required, if one wants to determine the 
polarization degree to an accuracy of order P. 

To evaluate the error A numerically, we assume that the Aj a are normally distributed, 
and that the projection of the angular momentum on the orbital plane ( relative to the 
rr-axis is uniformly distributed between and 2n. Figure 1 shows the distributions of A for 
10 6 random realizations in the case of a 1% polarization degree (P = 10~ 2 ). The SNR was 
assumed to be p = 100, 300, 500 or 1000, for which the fractions of the realizations with large 
errors (|A| > 0.5) are 80%, 44%, 20% and 1%, respectively. Therefore, p ~ 10P' 1 = 1000 is 
required to determine a 1% polarization degree. 

The LIGO interferometers are coaligned, and thus cannot by themselves determine the 
polarization degree of GWs, although a network consisting of LIGO and other interferom- 
eters could in principle do it. We assume here an interferometer identical to the LIGO 
detectors with an optimal orientation (rotated by tt/4 with respect to the LIGO interfer- 
ometers). If some fraction of GRB progenitors are double neutron star mergers, the SNR 
is given by p ave ~ 16(<i/100Mpc)~ 1 (e.g. eq. (16) in Kobayashi & Meszaros 2002) where 
d is the distance to the binary. This estimate was obtained from an average over differ- 
ent possible orientation of the source and interferometer. The SNR p is larger by factor 
of (5/4) [(1 + cos 6f + (1 - cos eyf' 2 . Then, in an optimal case, we can determine a 1% 
polarization degree up to d max ~ 7 Mpc. If we assume that ~ 1000 GRBs happen in a year 
within ~ 3000Mpc, the closest event in a year is, on average, at d ~ 300 Mpc. Since the 
position in the sky is random, the SNR is smaller by a factor of a/2/5 than in the direction 
of the z-axis. The LIGO is most sensitive around / ~ 150Hz and the sensitivity is about 
V foSh(fo) ~ 3 x 10~ 23 . Since d max is proportional to S h 1 ^ 2 & future detector with sensi- 
tivity ^/ foSh(fo) ~ 4 x 10 -25 would enable us to measure a polarization degree of 1% in a 
timescale of one year. 
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Frail et al. (2001) show that the observed distribution of jet opening angles (or viewing 
angles for the angle dependent jet model) is f b s oc 9j 3 ' 5 (9j > 0.05). Assuming this result, 
the probability that GWs from a GRB have a polarization degree larger than P = 10~ 2 
is given by ~ 3 x 10 -2 (P/10~ 2 )~ ' 4 . However, the f b s are estimated from observations of 
typical GRBs at redshifts z ~ 1, and one might fail to detect a significant fraction of the 
dim GRBs at large viewing angles. Therefore, when we study GWs from nearby GRBs, the 
probability that one observes GWs with a large polarization degree could be much larger 
than the above estimate. 



5. Discussion and Conclusions 

We have argued that GRB progenitors are likely to emit I = m = 2 gravitational 
waves (GWs), which are circularly polarized on the polar axis, while the + polarization 
dominates on the equatorial plane. Recent GRB studies suggest that the wide variation in 
the apparent luminosity of GRBs are caused by differences in the viewing angle, or possibly 
also in the jet opening angle. Since GRB jets are launched along the polar axis of GRB 
progenitors, correlations among the apparent luminosity of GRBs and the amplitude as well 
as the polarization degree of the GWs are expected. 

At a viewing angle larger than the jet opening angle 9j (which may be defined also 
in the case of a universal angle-dependent jet profile) the GRB 7-ray emission may not be 
detected. However, in such cases an "orphan" long-wavelength afterglow could be observed, 
which would be preceded by a pulse of GWs with a significant linearly polarized component. 
An expanding jet with an opening angle 9j behaves, as long as its Lorentz factor 7 > 9j l , 
as if it were part of a spherical shell, but relativistic beaming effect allows only observers at 
viewing angles < 9j to observe the emission from the jet. (By contrast, since GW emission 
is expected from the central engine itself, GWs are not subject to such extreme beaming). 
As the jet slows down and reaches 7 ~ 0J 1 , the jet begins to expand laterally, and its 
electromagnetic radiation begins to be observable over increasingly wider viewing angles. 
Since the opening angle increases as ~ 7" 1 oc t 1 / 2 (Sari et al 1999), at a viewing angle 
9 > 9j, the orphan afterglow begins to be observed (or peaks) at a time t p oc 9 2 after the 
detection of the GW burst. The polarization degree and the peak time should be correlated 
as P oc t 2 . 

A new type of fast transient source, called "X-ray flashes" , have recently been observed 
with the BeppoSAX satellite (Kippen et al. 2002). Apart from their large fraction of X-rays 
(~ 2 — 10 keV), the overall properties of these events are similar to those of GRBs. Recently, 
Yamazaki, Ioka and Nakamura (2002) and Woosley, Zhang & Heger (2002) suggested that 
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these events may be GRBs with large viewing angles. If this is the case, linearly polarized 
GWs should be observed prior to the X-ray flashes. The degree of polarization should be 
positively correlated with longer delays and with the softness of the X-ray flashes, which 
increase with angle. 

Since the degrees of linear and circular polarization depend on the viewing angle, a 
determination of the polarization degree would be a measure of the viewing angle. Such 
measurements, which are likely to require the advent of a future generation of detectors, 
could provide a new tool for estimating the absolute luminosity of GRBs, including its pho- 
ton component. By comparing the estimated absolute photon luminosity with the apparent 
luminosity, the distance to the source may be estimated independently of any redshift mea- 
surement. No optical afterglows have been found for about half of all the GRBs detected by 
BeppoSAX (the so called "dark GRBs"), and the present method would have the potential 
to help determine or constrain the distances to such dark GRBs. 

We thank L. S. Finn, B. Owen, I. Jones and P. Sutton for useful discussions. We 
acknowledge support through the Center for Gravitational Wave Physics, which is funded 
by NSF under cooperative agreement PHY 01-14375, and through NSF AST0098416 and 
NASA NAG5-9192. 
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Fig. I. — Distribution of errors A for a polarization degree P = 10 2 , signal to noise ratios 
p = 100 (thin solid), 300 (dashed), 500 (dashed dotted) and 1000 (thick solid), based on 10 6 
random realizations. The size of the bins is 0.1. 



